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Abstract 

Let {Xjk)j,k^i be an infinite array of i.i.d. complex random variables, with 
mean and variance 1. Let \ n i, . . . , X n n be the eigenvalues of (A=Xjk)\<j k<n- 
The strong circular law theorem states that with probability one, the empir- 
ical spectral distribution j L {S\ n 1 + • • • + $\ nn ) converges weakly as n — > oo 
to the uniform law over the unit disc {z E C; \z\ ^ 1}. In this short note, 
we provide an elementary argument that allows to add a deterministic matrix 
M to (X jk ) 1<jik<n provided that Tr (MM*) = 0(n 2 ) and rank(M) = 0(n a ) 
with a < 1. Conveniently, the argument is similar to the one used for the 
non-central version of Wigner's and Marchenko-Pastur theorems. 
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1 Introduction 

q 

For any square n x n matrix A with complex entries, let the complex eigenvalues 
Ai(A), . . . , A n (A) of A be labeled so that |Ai(A)| ^ • • • ^ |A„(A)|. The empirical 
spectral distribution of A is the discrete probability measure //a : = r Yl k =i ^ A fe( A )- 
We denote by Si(A) ^ • • • ^ s n (A) the singular values of A, i.e. the eigenvalues 
of the positive semi-definite Hermitian matrix \J A A* where A* is the conjugate- 
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transpose of A. The operator norm is si(A) = max^y = i ||Aa;|| 2 and the square 
Hilbert- Schmidt norm is || A|| 2 := s 1 {A) 2 + - • - + s n (A) 2 = Tr(AA*) = Y>lk=i \ A j,k\ 2 - 
Weyl's inequality |Ai(A)| 2 H h |A n (A)| 2 ^ sx(A) 2 H hs n (A) 2 ensures that the 
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second moment of ^tA is always bounded above by -||A|| . The following result was 
recently obtained by Tao and Vu [THl Corollary 1.15]. 

Theorem 1.1 (Circular law for central random matrices). Let (Xj k )j t k^i be i.i.d. 
complex random variables. Let {Mj k )j^\ be deterministic complex numbers. For 
every integer n^l, set X„ = (X jk )i^j,k^n and M n = (Mjk)^^. If 

• E[|X l5l | 2 ] = 1 andE[X 1A ] = 

• ||M n || 2 = 0(n 2 ) and rank(M n ) = 0(n a ) for some a < 1 
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then with probability one, Ai_^(x n +M n ) tends weakly as n — >■ oo to the uniform dis- 
tribution on the unit disc {z G C; \z\ 1} (known as the circular law). 

The aim of this note is to provide an alternative and elementary argument which 
reduces theorem ll.ll to the central case where M„ = for every n. Conveniently, the 
approach is close in spirit to the one used by Bai [3] for the derivation of Wigner's 
and Marchenko-Pastur theorems for non-central random matrices. 

This note was motivated by the study of random Markov matrices, including the 
Dirichlet Markov Ensemble [El E], for which a circular law theorem is conjectured. 
The initial version of this note was written before the apparition of [TH] , and provided 
for the first time a non-central version of the circular law theorem. The initial version 
was based on potential theoretic tools. For convenience, the present version makes 
use instead of the replacement principle borrowed from j!9j . 

Theorem 11.11 belongs to a sequence of works by many authors, including Mehta 
[13], Girko [TO], Silverstein [T2], Bai [2J, Edelman [S] Sniady [T7], Bai and Silverstein 
[I], Pan and Zhou [H], Gotze and Tikhomirov [11], and Tao and Vu |18j . 



Remark 1.2 (Constant case). Consider the case where the entries of M n are all 
equal to 1 in theorem \l.l[ We have then rank(M„) = 1 and Si(M„) = n. Suppose 
additionally that X lt i has finite fourth moment. Then, by Bai and Yin theorem |5]/, 
with probability one, lim^oo si(^X n ) = 2, and thus ^(X n + M n ) is a random 

bounded perturbation of the rank one symmetric matrix which has spectrum 

A„,(M n ) = • • • = A 2 (M n ) = and A x (M n ) = V™. 

From this observation, Silverstein fib!/ has shown, via perturbation techniques such 
as Bauer-Fike and Gerschgorin theorems, that with probability one, 



A,( -^(X n + M r 

n 



^2 + o(l) and 



A i ( ^(X n + M n ) | - s/'n 



< 2 + o(l). 



See also the work of Andrew J7]/. Also, with probability one, as n — > oo, the spectral 
radius |Ai(^(X n + M n ))| blows up while /-i^_(x n +M„) remains weakly localized. 



2 Reduction to the central case 

In order to show that theorem 11.11 reduces to the central case where M n = for 
every n ^ 1, it suffices to check the assumptions of the replacement principle of 
theorem 13. II with A n := ^X n and B n := ^(X n + M n ). By the strong law of large 
numbers and the assumption on ||M n ||, with probability one, 

-||A n || 2 + -||B n || 2 = 0(l). 

n n 

Next, by theorem (j3.4p and the first Borel-Cantelli lemma, for all z G C, with 
probability one, the random matrices A n — zl n and B n — zl n are invertible for large 
enough n. Let us define, for large enough n, the quantity 

1 1 

A n , 2 := - log |det (A„ - zln)] log |det (B„ - zln)]. 
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If we set ^ := V ^ {An _ zln){An _ zlnY and u n , z := ^ (BB _ aIn)(BB _ aIn) . then 

/•oo 

A n ^ = / log(t) d{ii n>z - v n>t )(t). 
Jo 

By the strong law of large numbers and the assumption on ||M n ||, for all z G C, 
with probability one, there exists a > such that 

max(si(A„ - zl n ), si(B n - zln)) < n a 

for large enough n. On the other hand, by theorem (13.41) and the first Borel-Cantelli 
lemma, for all z G C, with probability one, there exists b > such that 

min(s n (A n - zl n ), s„(B n - zl n )) ^ 

for large enough n. Therefore, with a n := n~ b and /3 n := n a , and large enough n, 

log(t) d(yU„ >2 - U n j)(t). 

n 

Let F n z and G n ^ z be the cumulative distribution functions of the real probability 
measures fj, HrZ and i/ HrZ . By lemma I3T31 and the assumption on rank(M ra ), for almost 
all z G C, with probability one, there exists e > such that 

IK* - G n J =0(n- £ ). 



Therefore, by lemma I372j we obtain, for almost all z G C, with probability one, 
|A n , z | < (log(A.) - log(a„))||F n , 2 - G n J\ = o(l). 



3 Tools 

This section gathers some tools used in our proof of theorem 11.11 By Green's theo- 
rem, for any complex polynomial P and smooth compactly supported / : C — >■ R, 

/ fdfi = — I A/log |P| dxdy 

where /i := 5a x + • — h <$a„ is the counting measure of the roots Ai, . . . , A n of P in C. 
Used for characteristic polynomials of random matrices, this identity provides, via 
dominated convergence arguments, the following theorem, see [T9| Theorem 2.1]. 

Theorem 3.1 (Replacement principle). Let (A n ) n >i and (B n ) n ^i be two sequences 
of complex random matrices where A n , B n are n x n, without any assumptions. If 

• with probability one ^||A n || 2 + ^||B n || 2 = 0(1) 

• for almost all z G C, with probability one, the random matrices A n — zl n and 
B n — zl n are invertible for large enough n 
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• for almost all z G C, with probability one, 

lim f - log [det (A n - zl n )\ - - log |det (B n - zl n 

then with probability one, ixx n — ^B n tends weakly to zero as n — > oo. 

The following lemma is a special case of the integration by parts formula for the 
Lebesgue-Stieltjes integral (with atoms). We give a short proof for convenience. 

Lemma 3.2 (Integration by parts). If ax, . . . , a n , bi, . . . , b n G [a, (3] C R, and 

and F v are the cumulative distribution functions of fi — -(5 ai + ••• + 8 an ) and 
v = + • • • + b~b„) respectively, then for any smooth f : [a, (3\ — > R, 

f(x)dfi(x)— / f{x)dv(x)= I f'(x)(F fM (x) — F v (x)) dx. 



In particular, when f is non decreasing, 

f(x)d fJ ,(x)- I f{x)dv{x) ^{m-fiaMF^-F^. 

J a 

Proof. One can assume by continuity that ai, . . . , a n ,bi, . . . ,b n are all different. We 
reorder oi, . . . , a n , 6 1} . . . , b n into c x ^ • • • ^ c 2n . For every 1 ^ k ^ 2n, set e k = +1 
if Cfc G {ai, . . . , a n } and e k = —1 if G . . . , b n }. We have 



pP pP Y n y 2n 

/ f(x)dfi(x)- /(x)di/(i) = -V(/(a i )-/(6i)) = -Ve fc /(c fc j 



By an Abel transform, we get by denoting S k = e\ + ■ — I- s k , 

2n 2n-l 

5^£fe/(c fc ) - - ^ S k (f(e k +i) - /(Cfc)) + S 2 nf(c 2 n)- 

k=l k=l 

Since F M — F u is constant and equal to S k on c k+ i[, 

S k (f(c k+1 ) - /(c fc )) = r +1 f'(x)(F,(x) - F u {x)) dx. 

It remains to notice that SWi = F^(c2 n ) — F u (c2 n ) = 0. □ 

The following lemma is a direct consequence of interlacing inequalities for singu- 
lar values obtained by Thompson J2U] in 1976. It was also obtained by Bai [5] 
and generalized by Benaych- Georges and Rao |6j. It is worthwhile to mention 
that it gives neither an upper bound for Si(B), . . . , s k (B) nor a lower bound for 
Sn-/t+i(B), . . . , s n (B) where k := rank(A — B), even in the case k = 1. 

Lemma 3.3 (Rank inequality). Let A and B be two n x m complex matrices. Let 
^/aa* > ^bb* ^ e the cumulative distribution functions of H^xa* an< ^ I^Vbb*- Then 

- ^VbbHL < ^ank(A - B). 
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The following theorem is due to Tao and Vu [HI Theorem 2.1], and is inspired 
from the work of Rudelson and Vershynin [15] . 

Theorem 3.4 (Polynomial bounds for smallest singular values). Let L be a probabil- 
ity distribution on C with finite and non-zero variance. For every constants A > 
and C\ > 0, there exists constants B > and C*2 > such that for every n x n 
random matrix X with i.i.d. entries of law L and every n x n deterministic matrix 
C with si(C) ^ n 1 , we have 

P(s n (X + C) ^ n- B ) ^ C 2 n- A . 
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